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Abstract 

Background: Investigating the host response in the early stage of influenza A virus (lAV) infection is of 
considerable interest. However, it is conceivable that effects due to the anesthesia and/or intranasal infection 
procedure nnight introduce artifacts. We therefore ainned to evaluate the effects of anesthesia and/or intranasal 
infection on transcription of selected pulmonary nnRNAs in two inbred nnouse strains with differential susceptibility 
to lAV infection. 

Results: DBA/2J and C57BLy6J nnice were evaluated in a tinne course experinnent in which lung tissue was sannpled 
after 6, 12, 18, 24, 48 and 120 h. After anesthesia with ketamine and xylazine, a suspension of mouse-adapted lAV strain 
PR8_l\/lun in 20 |jl sterile buffer, or 20 \a\ sterile buffer only, was instilled intranasally. The mice receiving anesthesia and 
PBS only were designated the "mock treatment" group. Pulmonary expression of 10 host mRNAs {Fos, Retnia, Irgl, 116, 
lllb, CxcllO, Statl, Ifng, Ifnl2, and Mxl) and viral hemagglutinin (HA) mRNA were determined at the designated time 
points. As expected, weight loss and viral replication were greater in the DBA/2J strain (which is more susceptible to 
lAV infection). Four mRNAs {Retnia, Irgl, 116, and CxcllO) were procedure-dependently regulated in DBA/2J mice 
between 6 and 24 h, and two (Retnia and 116) in C57BL/6J mice, although to a lesser extent. All 10 mRNAs rose 
after infection, but one (Fos) only in DBA/2J mice. These infection-dependent effects could be separated from 
procedure-dependent effects beginning around 12 h in DBA/2J and 18 h in C57BL/6J mice. The interferon-related 
mRNAs Statl, Ifng, Infl2, and Mxl were unaffected by mock treatment in either mouse strain. Mxl and Infl2 correlated 
best with HA mRNA expression (r = 0.97 and 0.93, respectively, in DBA/2J). 

Conclusions: These results demonstrate effects of the anesthesia and/or intranasal infection procedure on 
pulmonary gene expression, which are detectable between approximately 6 and 24 h post procedure and vary in 
intensity and temporal evolution depending on the mouse strain used. Mock infection controls should be included 
in all studies on pulmonary gene expression in the early phase of infection with lAV and, likely, other respiratory 
pathogens. 
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Background 

Due to ease of infection, animal rearing, and the avail- 
ability of genetically modified strains, using mouse 
models and viral strains adapted to the murine host has 
become an attractive approach to studying the mamma- 
lian response to influenza A virus (lAV) infection. Re- 
cently, a substantial amount of information has been 
obtained regarding gene expression changes at various 
stages of infection in this model [1-3]. These authors 
showed that the genetic background of different mouse 
strains strongly influences the susceptibility to lAV. For 
instance, DBA/2J mice were highly susceptibility to the 
PR8_Mun lAV (HlNl) strain, whereas C57BL/6J mice 
were more resistant, as evidenced by less weight loss 
and an approx. 1000-fold higher viral LD50. Conversely, 
viral load was significantly higher in the DBA/2J strain, 
which also mounted a hyper-inflammatory response with 
much stronger up-regulation of many immune response- 
dependent genes. As exemplified by the aforementioned 
studies, most work in murine models of lAV infection has 
focused on time points during or after established infec- 
tion (1 day up to 60 days), and very little attention has 
been paid to the first 24 hours (h). Nevertheless, critical 
aspects of the host response to early steps in viral attach- 
ment and entry could conceivably be studied during this 
early time window. However, due to the temporal proxim- 
ity to the technical and pharmacological manipulations 
surrounding the infection process, it is conceivable that 
both the administration of the anesthetic and the physical 
and physiological stress from intranasal installation of the 
inoculate would lead to artifactual signals that are unre- 
lated to the virus-host interaction. We have therefore ana- 
lyzed changes in pulmonary gene expression in a 5 -day 
time course featuring frequent measurements in the first 
24 h, comparing results obtained from mice infected with 
lAV or exposed to vehicle only ("mock infection"). We 
find effects on pulmonary gene expression that can be 
clearly ascribed to the anesthesia/infection procedure, 
which are detectable as early as 6 h post treatment and 
differ between the two mouse strains in terms of magni- 
tude and temporal evolution. 

Methods 

Sample preparation 

Female 12- 13- week- old C57BL/6J and DBA/2J mice (n = 
5-8 per time point and treatment) and mouse-adapted 
lAV strain variant PR8_Mun (Institute of Molecular Vir- 
ology, University of Muenster, Germany), which is closely 
related to A/Puerto Rico/8/34, were used. Mice were 
weighed on day 0 just before induction of anesthesia and 
on each subsequent day. Infections were essentially carried 
out as described previously [1]. Briefly, mice were anesthe- 
tized by intra-peritoneal injection of 10 [A per g body 
weight of a stock solution of 0.5 ml ketamine (50 mg/ml. 



Invesa Arzneimittel GmbH, Freiburg, Germany), 0.5 ml 
2% xylazine hydrochloride (Bayer Health-Care, Leverku- 
sen, Germany) and 9 ml sterile NaCl 0.9% (Delta-Select 
GmbH, Dreieich, Germany). For intranasal infection, a 
viral dose of 2 x 10^ focus forming units (ffu) of PR8_Mun 
(propagated in embryonated chicken eggs) was adminis- 
tered in a total volume of 20 [A sterile phosphate-buffered 
saline (PBS). During the infection procedure, mice were 
held in the upright position and additional anesthetic was 
reinjected as needed. Mock treatment was identical to real 
anesthesia/infections except that vehicle only (sterile PBS), 
not containing virus, was used for intranasal instillation. 
Mice were killed by CO2 asphyxiation at 6, 12, 18, 24, 48, 
and 120 h with respect to infection or mock treatment. 
Untreated mice were used as t = 0 h control. Lungs were 
removed from the cadavers by cutting the main bronchus 
and were washed in RNAlater RNA Stabilization Reagent 
(Qiagen Inc, Valencia, CA, USA) immediately after re- 
moval. After transfer into a new tube containing 2 ml 
RNAlater, lungs were stored overnight at 4°C and then at 
-20°C until further use. All animal work was approved by 
an external committee according to the regulations on 
animal welfare of the Federal Republic of Germany. 

RNA isolation and qRT-PCR 

Lungs were homogenized in 4 ml RLT buffer (Qiagen) 
containing 40 [A p-mercaptoethanol and stored at -80°C 
in 450 [A aliquots. After thawing, 450 [A of this suspen- 
sion was mixed with 700 [A Qiazol (Qiagen), and all fur- 
ther steps of total RNA isolation were performed with 
the miRNeasy kit (Qiagen) according to the manufac- 
turer s recommendations. Real-time RT-PCR (qRT-PCR) 
was performed with a LightCycler 480 (La Roche AG, 
Basel, Switzerland) in 96 well plates in 20 (il reaction 
volumes, using 15 ng cDNA (miScript Reverse Tran- 
scription Kit, QuantiTect SYBR Green PGR Kit) and 
primers specific for the following targets: the immediate 
early gene FBJ osteoscarcoma oncogene (Fos), resistin 
like a {Retnla), immune-responsive gene 1 {Irgl), interleu- 
kin 6 {116), interleukin 1|3 {Illb), the chemokine (C-X-C 
motif) ligand 10 {CxcllO), four genes related to interferon 
pathways (the transcription factor signal transducer and 
activator of transcription 1 {Statl), interferon y {Ifng), 
interferon \2 {Ifnl2, aka Il28a), and myxovirus (influenza 
virus) resistance 1 {Mxl)X and lAV hemagglutinin (HA), 
Quantitect Primer Assays (Qiagen) were used for all tar- 
gets except Ifhl2 and HA. Primers for amplification of 
Ifhl2 were designed using exon-spanning regions of the 
NCBI [4] sequence (Tanta_Mus_Ifnl2-F: 5ctgcttgagaagg 
acctgagg'3, Tanta_Mus_Ifnl2-R: 5ctcagtgtatgaagaggctggc'3). 
Primer sequences for HA mRNA amplification were pub- 
lished previously [3]. Mouse Genome Informatics (MGI) 
gene symbols and names were used for all genes [5]. The 
arithmetic mean of the Ct values of |3 actin {Actb) and 
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ribosomal protein L4 {Rpl4) was used as internal reference 
for normalization. 

Data analysis 

Data were analyzed using the R environment and pro- 
gramming code [6]. qRT-PCR data points with Q >40, 
corresponding to lack of detection of a target due to tech- 
nical failure or lack of expression, were assigned a Q of 
40. We removed technical outliers in AQ values using the 
maximum normed residual test (Grubbs' test) to detect 
outliers for each condition with a threshold of p <0.05. A 
median of 5 (range, 3-8) biological replicates were avail- 
able for each data point after outlier removal. ANOVA 
was used for testing of trends throughout time series, 
adjusting p values for false discovery rate (FDR). For pair- 
wise comparisons, we used Tukeys Honest Significant 



Differences Test for homogeneous variances and Dunnett s 
Modified Tukey-Kramer Pairwise Multiple Comparison 
Test for heterogeneous variances (Levene s test for variance 
equality). We used a significance threshold of p <0.05. 

Results 

Changes in body weight 

Weight loss in infected DBA/2J mice became manifest 
by day 2 (Figure lA) and reached approx. 25% by day 5. 
As shown previously [3], weight loss in the infected 
C57BL/6J mice was less pronounced, as reflected in a 
mere 5 - 10% reduction by day 4 - 5. In both strains, 
statistically significant differences between infected and 
mock treated mice were observed by day 3. Mock- 
treated mice showed no significant weight loss at any 
time point. Thus, there was no significant effect of the 
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Figure 1 Weight loss and expression of lAV HA mRNA throughout the 5-day time course after mock treatment or infection with lAV 
strain PR8_MUN as outlined in the Methods section. A. Weight loss, expressed as the percentage of body weight measured at t = 0 h before 
administration of anesthesia. No mice had to be l<illed because of >30% weight loss. B. Relative quantification of lAV HA mRNA in mouse lung by 
qRT-PCR in the 5-day time course shown in panel A. dCt refers to Ctreference - Cttarget mRNA/ where Ctreference corresponds to the arithmetic mean 
of the Ct values of Actb and Rpl4. Solid lines, infection; interrupted lines, mock treatment. Left panels, DBA/2J strain; right panels, C57BI76J strain. 
Note that the x-axes of the panels are based on different scales. ^ p <0.05 for difference with respect to t = 0 h; i, p <0.05 for difference between 
mock-treated and infected mice at the given time point (Jukey's test). 
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anesthesia/infection procedure on body weight in either 
mouse strain. 



Viral replication 

qRT-PCR revealed a brisk rise of mRNA encoding lAV HA 
in lungs of both mouse strains after infection (Figure IB). 
HA mRNA was detected at low levels as early as 6 h in 
both strains, followed by a rapid rise that peaked at 48 h 
and 120 h in DBA/2J and C57BL/6J mice, respectively. HA 
mRNA levels were significantly higher in DBA/2J than in 
C57BL/6J starting around 12 h. As expected, HA mRNA 
was not detected in the mock treated mice. 



Principle component analysis of pulmonary expression of 
host-encoded mRNAs 

A cluster containing infected and mock treated time 
points could be identified easily in both mouse strains 
(Figure 2). A separation between infected and mock- 
treated samples became evident at 18 h in both mouse 
strains, as indicated by the lines in Figure 2. Marked 
step-offs between 24 and 48 h were seen in both strains. 
Consistent with the continuing rise of HA mRNA in the 
C57BL/B6 strain between 48 and 120 h the 120 h time 
point localized beyond the 48 h time point. In contrast, 
in the DBA/2J strain HA mRNA declined between 48 
and 120 h, and the 120 h time point localized between 
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Figure 2 Principal component analysis of pulmonary expression of the 10 host-encoded mRNAs in mock-treated and infected DBA/2J 
and C57BL/6J mice in the 5-day time course of lAV infection. mRNA levels of Fos, Retnlo, Irgl, 116, lllb, CxcllO, Stotl, Ifng, Ifnl2, and Mxl were 
determined by qRT-PCR as outlined in the Methods section, using Actb and Rpl4 mRNA expression for internal normalization. Each dot refers to 
the mean value of mice of one treatment as outlined in the legend adjacent to the box. The number inside each dot identifies the time (h) 
elapsed since t = 0 h. A. Results obtained with the DBA/2J strain. B. Results obtained with the C57BL/6J strain. Including the third component did 
not lead to further discrimination (data not shown). 
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24 and 48 h in the PCA plot. In both strains, the t = 
48 h and 120 h mock treated mice locaUzed far away 
from the infected t = 48 and 120 h mice. Thus, this ana- 
lysis showed that in both mouse strains procedure- 
related effects did not contribute significantly to the 
gene expression changes observed in infection beyond 
24 h and that the evolution of gene expression reprogram- 
ming correlated well with transcription of HA mRNA in 
both strains. 

Pulmonary expression of individual host-encoded mRNAs 

Results are shown in Figure 3. All 10 host mRNAs exhib- 
ited at least some evidence of regulation throughout the 
time course (ANOVA). Four mRNAs were also signifi- 
cantly regulated in response to mock treatment, but two of 
these {CxcllO and Irgl) were regulated only in the DBA/2 J 
strain. Fos, Illb, Statl, Ifhg, Ifnl2, and Mxl mRNAs were 
not regulated by mock treatment in either strain. 

Fos mRNA was expressed at low level and increased 
significantly after 48 h in the infected DBA/2J mice only 
(panel A). There was no evidence for mock treatment- 
dependent regulation of this mRNA in either mouse 
strain (panels A and B). The apparent tendency of an in- 
crease in infected C57BL/6J mice toward the later time 
points was not significant. Retnla mRNA increased sig- 
nificantly at all time points in infected DBA/2J mice. 
ANOVA revealed a significant effect of the mock treat- 
ment, with an effect size resembling that of infection be- 
tween 6 and 24 h, suggesting procedure-dependent 
regulation through 24 h and infection-dependent regula- 
tion thereafter (panel C). In infected C57BL/6J mice, 
Retnla increased between 18 h and 120 h (panel D). In 
mock treated mice of this strain, expression increased 
steadily at the early time points, with significant regula- 
tion at 18 h and 24 h. This suggested a procedure- 
dependent regulation of Retnla resembling that observed 
in the DBA/2J mice. Irgl mRNA expression increased in 
mock-treated DBA/2 J mice between 6 and 18 h (panel 
E). This gene was up-regulated in DBA/2 J infected mice 
at all time points, reaching a maximal 630-fold induction 
on day 2. In the C57BL/6J strain there was no increase 
in Irgl due to mock treatment, and the infection- 
dependent increase was less pronounced, reaching a 
max. 150-fold induction at 120 h. 116 mRNA increased 
in both strains beginning 6 h after infection or mock 
treatment, with stronger regulation being observed in 
the DBA/2 J mice (panel G). In DBA/2 J mice the mock 
treatment effect declined towards 18 h, and clear differ- 
ences between infected and mock treated mice became 
apparent at 24 h. In the C57BL/6J mice, an infection- 
dependent rise in 116 mRNA was observed somewhat 
later (t = 48 h) (panel H). II lb mRNA increased in in- 
fected mice of both strains at 48 h and 120 h, and there 
was a tendency (p at 6 h = 0.09) toward a mock 



treatment effect between 6 and 18 h in the DBA/2 J strain 
(panel I). CxcllO mRNA was up-regulated in DBA/2J 
mock-treated mice at 6 h (panel K), whereas it was not af- 
fected by mock treatment in the C57BL/6J mice (panel L). 
In both mouse strains CxcllO mRNA was significantly ele- 
vated in the infected mice, beginning at 6 h in the DBA/2J 
and at 18 h in C57BL/6J. Statl expression was not affected 
by mock treatment in DBA/2J mice, but there was a slight 
trend (statistically not significant) for up-regulation in 
C57BL/6J mice. An infection-dependent up-regulation be- 
came apparent at 24 h and 48 h in DBA/2J and C57BL/6J 
mice, respectively. Similar to Static Ifng was up-regulated 
in both mouse strains beginning around 48 h, and there 
was no evidence for regulation due to the infection pro- 
cedure. Ifnl2 was not detected (Q > 40) in about 40% of 
untreated and mock treated DBA/2J mice; fold change 
values therefore represent an underestimation (panel Q). 
A significant rise after infection became apparent at 48 h, 
reaching a mean Q of 26.3. In C57BL/6J mice, it was not 
detected in about 80% of the untreated and mock treated 
samples, suggesting a lower baseline expression than in 
DBA/2 J (panel R). A first significant infection-dependent 
regulation was observed at 18 h, where Ifnl2 was detected 
in all DBA/2J and four of five C57BL/6J samples. Ifnl2 was 
detected in all 24 h samples (Q = approx. 33) and contin- 
ued to rise through 120 h. There was no evidence for a 
mock treatment effect on Ifnl2 in either mouse strain. 
Mxl mRNA expression (panels S and T) was not regu- 
lated in response to mock treatment in either strain. In 
contrast, it strongly increased during infection, showing a 
significant rise beginning at 12 h in DBA/ 2 J and 24 h in 
C57BL/6J mice. 

Taken together, these results allow classifying the ana- 
lyzed genes into three groups: (1) genes that were regu- 
lated in response to mock treatment and infection in 
both strains {Retnla, 116), (2) genes that were regulated 
in response to both mock treatment and infection in the 
DBA/2 J strain only {Irgl, CxcllO), and (3) those whose ex- 
pression changed in response to infection only {Fos, Illb, 
Statl, Ifng, Ifnl2, and Mxl), Of note, the latter group con- 
tained all four interferon pathway-related mRNAs. 

Correlation with lAV HA mRNA 

Expression of the 10 host mRNAs was then correlated 
with HA mRNA expression (Table 1). Overall, correla- 
tions were higher in the DBA/2 J strain. Only Illb corre- 
lated more strongly in C57BL/6J than in DBA/2J. Mxl 
and Ifnl2 mRNA levels correlated best with HA mRNA 
expression in both strains, whereas Fos mRNA was the 
only one that did not correlate with HA mRNA. 

Regulation across all 10 target mRNAs 

Results are summarized in Figure 4. Considering regula- 
tion across all 10 target mRNAs combined, we detected 
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(See figure on previous page.) 

Figure 3 Expression chianges in moclc-treated and infected DBA/2J (left column) and C57BL/6J (right column) mice in the 10 target host 
mRNAs in the 5-day time course of lAV infection. Analysis of tine same data set as used for Figure 2. Panels show fold change expression, 
determined by qRT-PCR, of Fos (A, B), Retnio (C D), Irgl {E, F), 116 [G, H), lllb {\, J), CxclW {K, L), Stotl {N\, N), Ifng {O, P), Ifnl2 [Q, R), and Mxl 
{S, T) mRNA. Fold change data of Ifnl2 represent an underestimation, as a Ct of 40 was assigned to all samples with Ct >40 (see Methods section). 
Solid lines, mice intranasally infected with 1x10^ ffu of lAV strain PR8_Mun; interrupted lines, mice undergoing the same anesthesia/infection 
procedure except that buffer only, not containing virus, was used for intranasal installation (mock treatment). *, p <0.05 for difference between 
infected mice at the given time point with respect to t = 0 h; A, p <0.05 for difference between mock treated and control (t = 0 h) mice; 
p <0.05 for difference between infected and mock treated mice at the given time point. All p values were determined with Tukey's test. 



a significant up-regulation at all time points after 0 h in 
infected DBA/2J mice (Dunnetts Modified Tukey- 
Kramer Pairwise Multiple Comparison Test). Among 
mock treated DBA/2J mice, an up-regulation was ob- 
served at 6, 18 and 24 h post treatment. The strongest 
effect was detected at 6 h (mean fold increase, 2.9; CI = 
1.6-5.4) which nearly equaled the regulation in infected 
mice (mean fold increase, 2.7; CI = 1.5-4.7). A significant 
difference between infected and mock-treated DBA/2J 
mice could be discerned by ANOVA beginning at 12 h, 
but a contribution of a procedure-related effect to 
mRNA expression in the infected mice could be ex- 
cluded only from 48 h onward. Messenger RNA up- 
regulation peaked at 48 h and began to decline by 120 h. 
In the C57BL/6J strain, overall up-regulation was less 
than in the DBA/2J strain. In this strain, the expression 
change at 6 h seemed to be due to the anesthesia/infec- 
tion procedure in both infected and mock-treated mice, 
as fold induction was nearly identical in both (mean fold 
induction, 1.6; Cljnf = 0.98-2.6 and CiMock = 0.84-2.9). 
As in the DBA/2J strain, a procedure-dependent effect 
seemed to persist through 24 h (CiMock = 0.97-2.23). 
Infection-dependent mRNA up-regulation first became 
manifest at 18 h and continued to rise between 48 and 
120 h. Thus, this analysis supported the kinetics of evo- 
lution of infection-dependent transcription in both 
mouse strains, as revealed by the PC A (Figure 2). In 
addition, consistent with the results shown in Figure 3, 
it showed that procedure-dependent effects occurred be- 
fore 48 h and were more pronounced in the DBA/2J 
strain. 

Discussion 

This analysis of sequential changes in pulmonary expres- 
sion of several mRNAs after real or simulated lAV infec- 
tion revealed effects that can be ascribed to anesthesia 
and/or the intranasal inoculation procedure. The results 
clearly demonstrate that the appropriate control group 
treated with a simulated anesthesia/infection should al- 
ways be included in studies of lAV infection in mice that 
cover approximately the first 24 h post infection. 

What might be the underlying pathophysiological 
mechanisms? Anesthesia is known to influence cytokine 
expression in humans, but actually appears to have an 



anti-inflammatory effect as, for instance, suggested by 
reduction of circulating 116 levels [7-9]. The intranasal 
infection procedure appears to be a more likely candi- 
date. Despite the relatively small volume of 20 [A that is 
used and the near physiological properties of PBS, we 
consider it likely that entry of PBS into the airway cre- 
ates a stress response similar to that observed after fluid 
aspiration, including at least focal pulmonary hypoxia 
due to bronchospasm. Responsible mechanisms may 
both relate to stimulation of nerve endings in the airway 
epithelium and direct noxious stimulation of airway epi- 
thelial cells. Indeed, except for Irgl, three of the four 
mRNAs whose expression was regulated in response to 
mock treatment are known to be induced during a stress 
response or hypoxia at the cellular or tissue level 
{Retnla: [10]; 116: e.g., [11]; CxcllO: [12]). The fourth 
one, Irgl, is preferentially expressed in macrophages, is 
strongly induced during macrophage activation, and lo- 
calizes to mitochondria [13,14]. Its expression in stress 
or hypoxia has not been examined, and it would there- 
fore be interesting to test whether it plays a role in these 
processes. The four interferon related genes {Statl, Ifng, 
Ifnl2 and Mxl) were clearly induced in infected mice 
only. They have all been shown to play important roles 
in the host defense against lAV and other viruses 



Table 1 Correlations of pulmonary expression of 1 0 target 
mRNAs with HA mRNA^ 



mRNA 


DBA/2J 


C57BL/6J 


Mxl 


0.97*** 


0.89*** 


Ifnl2 


0.93*** 


0.87*** 


CxcllO 


0.92*** 


0.87*** 


Stotl 


0.90*** 


0.86*** 


lie 


0.80*** 


0.68*** 


Ifng 


0.70** 


0.62** 


Irgl 


0.76*** 


0.72*** 


Retnla 


0.62** 


0.63*** 


lllb 


0.53* 


0.71*** 


Fos 


0.39 


0.16 



Values correspond to Spearman correlation coefficient in mouse strains 
infected with lAV, sorted by decreasing values in DBA/2J mice. P values (FDR 
adjusted): <0.001; **, <0.01; *, <0.05. 
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Figure 4 Overall mean fold changes in mRNA expression throughout the time course. A. Mean expression changes in mocl<-treated and 
infected DBA/2J and C57BI76J mice across all 10 target host mRNAs in the 5-day time course of lAV infection. The analysis is based on the same 
data set as used for Figures 2 and 3. Mean fold change values and 95% confidence intervals (vertical lines) were calculated with the Dunnett's 
Modified Tukey-Kramer test, using the dCt values (qRT-PCR) of all 10 host-encoded mRNAs as input. B. Schematic representation of the results 
shown in panel A. As reflected by the thickness of the lines, overall changes are more pronounced in the DBA/2J strain. Procedure-dependent 
effects are evident between 6 and 24 h in both strains, but infection-related changes begin to evolve and peak earlier in the DBA/2J strain. 



[15-17]. The Mxl gene is nonfunctional (truncated) in 
certain mouse strains including DBA/2J and C57BL/6J, 
but even the nonfunctional murine form is fully inter- 
feron inducible [18], suggesting that it does reflect the 
anti-influenza interferon response of the DBA/2J and 
C57BL/6J mice. Among these four genes, only Statl has 
been shown to be regulated by stress or hypoxia [19,20]. 
Interestingly, it was not affected by the mock treatment 
in the presented study, perhaps because its sensitivity to 
regulation in this murine model is not high enough to 
respond to any stress/hypoxia due to the mock treat- 
ment. Indeed, its upregulation in the infection was much 
smaller compared to the other three interferon-related 
genes. Thus, the observation that expression of these 
four interferon-related mRNAs was not affected by the 
mock treatment supports the aforementioned notion 
that the procedure-associated effects in this model relate 
to a stress response that can be functionally separated 
from the antiviral response. 

Differences between the two mouse strains 

Differences were observed in the magnitude of the re- 
sponse to both mock treatment and viral infection. The 



fact that both procedure and infection-related responses 
were more vigorous in the DBA/2J mice agrees with the 
previously described overall stronger inflammation in 
this strain during lAV infection [1]. This may reflect a 
greater intrinsic propensity to inflammation, but also the 
higher rate of viral replication in this strain. We favor a 
combination of both models, as the procedure-dependent 
effects, too, were brisker in the DBA/2J mice. 

Limitations 

The relatively small sample size represents a limitation 
of this study. Nonetheless, statistical significance was 
reached for several variables. A larger sample size would 
likely reveal additional significant changes, such as 
procedure-dependent regulation of Illb, at least in the 
DBA/2J strain, in which there currently is a tendency to- 
ward significance (mean fold increase at 6 h in mock- 
treated mice = 2.8; p = 0.09). In addition, the small num- 
ber of target mRNAs does not represent overall gene 
expression in the lung. Other methods such as RNA 
deep sequencing would likely reveal genes showing an 
earlier response to lAV infection or a longer persistence 
of procedure-dependent effects. 
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Conclusions 

Despite the aforementioned limitations, the presented 
results clearly show that the manipulations surrounding 
the infection procedure can affect pulmonary gene ex- 
pression in a host strain-dependent manner for approx. 
24 h. Thus, "mock treatment" controls should be in- 
cluded in all murine studies on lAV infection where 
measurements are to be taken within approx. the first 
24 h. Likewise, such controls are likely needed in similar 
studies with other viral and non-viral respiratory patho- 
gens. The degree to which a mock effect might interfere 
with infection-dependent data would depend on the kin- 
etics of replication and spread of the respective infection 
in the host animal. 
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